CALL OPTION PRICES BASED ON BESSEL PROCESSES 



JU-YI YENW-t 2 ) AND MARC YOR( 3 )'W 



Abstract. As a complement to some recent work by Pal and Protter [9], 
we show that the call option prices associated with the Bessel strict local 
martingales are integrable over time, and we discuss the probability densities 
obtained thus. 



1. Introduction: some general remarks 

1.1. Let (M t ,t > 0) denote a continuous local martingale, taking values in R + . To 
any K > 0, we associate the process (K — M t ) + , t > 0. It is not difficult to show, 
after localizing M, that this process (K — M t ) + is a (bounded) submartingale, and, 
as a consequence, the function: 



(t) = E[(K-M t ) + ], t>0 



is increasing, and bounded (by K). The study of such functions, considered (essen- 
tially) as distribution functions, has been the subject of the Bachelier Course [HE], 
given by the second author. In particular, if M t — ► 0, there is the formula 



L- 



E [F t (K - M t )+\ = KE [Ftl (gK < t) ] (1) 

which is valid for every F t > 0, (Tt) measurable, and Gk — sup{£ : M t — K}. See 
also Madan-Roynette-Yor [H H [H]. 

1.2. The present paper is devoted to the study of the functions: 

m { K \t) = E[(M t -K) + ] =E[(K-M t y], t > 0, 

which play an important role in option pricing, as (wiL- (t)) is the European 
call price with strike K, and maturity t, associated with the local martingale 
(Mt). If (M t ) is a "true" martingale, then {(Mf — K ) + } is a submartingale, hence 
(rn- K \t),t > 0) is increasing. On the other hand, if {M t ,t > 0) is a strict local 
martingale, that is: a local martingale, which is not a martingale, then the function 
(m' K •*) is not in general increasing, or even monotone. 

1.3. The most well-known example of a strict local martingale is M t = 1/Rt, 
where (Rt,t > 0) denotes the BES(3) process, starting from 1, or, by scaling, 
equivalently from any r > 0. Then, the study of (m^\t)) in this particular case 
has been undertaken in a remarkable paper by S. Pal and P. Protter 0; the results 
of which have strongly motivated the present paper. 
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In the present paper, we take up again the study of this function {rrS K \t)) in 
this particular case; we show that: 

dt m K ^ (t) < co. 

Hence, up to a multiplicative constant (m^ K (t), t > 0) is a probability density on 
R+ ; we identify the Laplace transform of this probability, and describe it as the law 
of a certain random variable defined uniquely in terms of BES(3) process. This is 
done thanks to the Doob /i-transform understanding of BES(3) (from Brownian 
motion, killed when hitting 0), combined with general identity JT]). We refer the 
reader to Section [2] for precise statements. In Section [3j we develop the same kind 
of study but this time with M t = l/i?^ , t > 0, where (Rt, t > 0) denotes 
the BES(5) process, starting from 1. In Section [4l we present the graphs of the 
corresponding functions {rrS K (t), t > 0). 

1.4. To summarize, the main point of this work is to use the interpretation of the 
generalized Black-Scholes quantities in terms of last passage times (formula |T])) in 
the framework of Bessel processes in order to derive fine properties of the call option 
process, as a function of maturity, written for the strict local Bessel martingales. 

2. Some results about r (3) {t) = e[ 3) [(-^ - K) + ], t > 

2.1. In this section, we change notation slightly: (X t , t > 0) denotes the canonical 
process on C(R+, R+), W x is Wiener measure such that W X (X = x) = 1, and Pi 
is the law of the BES(3) process starting from x. In fact, we shall only consider 
x = 1 (except mentioned otherwise). 

(31 

2.2. Here are our 3 main results concerning the functions r K (t). 
Proposition 1. The following holds: 

(i) ( t - co): rf (t) ~ yft; rf(t) ~ 
(ii) (t - 0): rf(t) - (1 - K)+; rf\t) ~ 7f . 



An important consequence of Proposition [I] (i), is that, for K > 0, the function 
{r K {t),t > 0) is integrable over R + , hence it is, up to a multiplicative constant, a 
density of probability on R + . We now describe this probability. 

Proposition 2. 

(i) The function {^K 2 r%\t),t> 0) is a probability density on R+ . 

(ii) It is the density of 

A$^(9i-T k )+T k U (2) 

where, on the RHS of $Sj), k — l/K, the variables g\, Tfc, and U are 
independent, g\ = sup{i : X t — 1}, T k is the size-biased sampling of 
T k = inf{t : X t = fcfl with gi and T k defined with respect to Pq , and 
finally U is uniform on [0, 1]. 

As a further description of the law of A^- , we present its Laplace transform. 



: That is: T k satisfies E[f(T k )} = 3K 2 E[f(T k )T k ] for every / : K + -+ K+, Borel. 
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Proposition 3. The Laplace transform of A K is: (we use again k = 1/K) 

'sinh(\/2Afc) \ 



E 



2Xk 



2.3. Proofs of Propositions 0} H El 



2.3.1. The main ingredients of these proofs are the following: 

(i) the Doob /i-process relationship between Brownian motion and BES(3), 
which may be written as: 

(ii) the particular instance of formula fl} with M t = X t AT a , under W\\ 

Wi (F t (k - X tATo ) + ) = kWi (F t l hk < t )) 

with 74, = sup{i < T : X t = k}. 

(iii) the time reversal result: {X To ^ tl t < T ) under W% is distributed as 
[X tl t<gi) under P (3) . 

2.3.2. Thanks to the preceding points, we may now obtain interesting description 

(3) 

of r K (t) in terms of first and last passage times. In fact, we obtain: 

rf{t)=W 1 { lk <t<T ) (3) 
or, equivalently, from the time reversal result in (iii): 

r${t) = Pl?\g l >t)-P£\g 1 >t). (4) 

Proof of ©: 

Combining (i) and (ii) above, we obtain: 

rfit) = W 1 ({l-KX t ) + l it<To) ) 

= ^((t-IwJ+l,^)), withfc = i 

= KkWx (Tfc < t < T ) 

= Wi (7fe < t < T ) , which is ©. 

2.3.3. Proof of Proposition [2] 

a) We deduce from ([3|) that: 

dt rP(t) = Wi (T - 7it ) = E^ 3) (T k ), 

by time reversal. Using the fact that (R 2 — 3t, i > 0) is a Pq 3 ^ -martingale, 
we get: 

3^ 3) (T fe ) = fe 2 . 

Hence, the constant c we were seeking is: c = p-, and (t) = 3K 2 r^p (t) 

is a probability density on M + . 
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b) In order to identify a random variable Aj^ with distribution 3K 2 r^ (i), we 
go back to ^ and we get, for any / : R + — > R + , Borel: 



/ [To 

dt f(t) 3K 2 r^\t) = 3K 2 W! ( / /(i) 

ik 



= 3K 2 W X [f ( 7fe + (T - Tfc) (7) (T - 7fc )] 

= 3X 2 T< 3) [/ (( 5l - T fe ) + T k U) T fe ] , 

where U is uniform on [0, 1], independent of Tt and g\. Hence using the 
notation of T k for the size-biased sampling of T k , we get that: 



E 



(/(Ak 3) 



= E, 



(3) 



f([9i-T k ) +T k U 



□ 



2.3.4. Proof of Proposition [3 

From formula (|3|) again, this Laplace transform is: 



3K 2 W l (^J dte 

The result now follows from: 

Wi (e~ AT °) = e" 

2.3.5. Proof of Proposition [T] 

(i) As t — > oo, we have: 
J3), 



-A/ 



. 2 A 



Wi (e 



2Xk 



sinh(V2Afc) 



□ 



,(3) m _ p (3) / ± \ _ ± f 1 
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The equivalent, for K > 0, of r^- (t), as i — > oo, is simply the particular 
case: 6 = 3 of the result in (i) of Proposition [4j 
(ii) The first statement follows from the convergence in L 1 of -j- to 1 as £ — > 0; 
For the second statement, we use: 

r[ 3) (t) = W x ((1 - X t ) + l (t<To) ) 

= ((1 - X t ) + ) - Wi ((1 - X t ) + l(T <t)) ■ 

The first term is: VF (A t + ) = 
jorized by 



j-, whereas the second one may be ma- 



W 1 (T <t) - W (T X <t) = W (\Xi\ > 

4 



(Note that, in fact, it is not necessary to use Proposition!]] to obtain Propositions^ 
and [31 however, it is the result of Proposition [TJ (i), which led us to study At?). 

2.3.6. A more direct proof of the estimate: r^'(t) = O (7572), as t — > 00. 

Since this estimate plays quite some role in our paper, it seems of some interest to 

look for a simple proof of it. We note that: 



-K 



(3) 



- K 



since, from the additivity property of squares of Bessel processes (Shiga- Watanabe 
|14j). a Bessel process with dimension S, starting from a > 0, dominates stochasti- 
cally a Bessel process with dimension 5, starting from 0. Now we have: 



E, 



(3) 



1 

Rt 



K 



E 



(3) 



( 1 X + 

- K 



\ViRi 

1/VtK 



(by scaling) 



dr r 2 e- r I 2 



< 





c 

VtJo 

c 



K 



for some universal constant c 



1/VtK 



dr 



-r 2 /2 



Vt 



1 — exp 



2tK 2 



t-oo 2t 3 / 2 K 2 ' 



3. Extending the previous results to: r^(t) = [(-^=2 - K) + ] , 5 > 2 



3.1. In this section, (X tl t > 0) still denotes the canonical process on 0^+,^+ ), 
and we consider P^ S) the law of the BES{8) process, starting from x (which, again, 
will be taken mainly equal to 1). 

3.2. As a parallel to Section ® we offer 3 results concerning the function r^) . We 
note v = I — 1 



-.(«) 



Proposition 4. The following holds: 

(1) (t -> 00): r^ 5 (t) ~ r 2 "T{i+v) > r x ) (*) ~ FTT' mth C )k = 2 "+i( I /+i)r(i/)A 

(ii) (t - 0): r£>(t) - (1 - K)+; r[ S) (t) ~ V^)^ 

Proposition 5. 

(i) TTie function is a probability density on R + . 



(ii) It is the density of 



A<? ( = W) ( 9 i-T fe )+f fe [/ 



(5) 



where, on the RHS of fc = l/K 1 /^ 2 \ the variables g-y, Tk, and U 
are independent, g\ = sup{t : X t = 1}, Tk is the size-biased sampling of 



Tfe = inf{£ : X t — k}, with g\ and Tf. defined with respect to P , and 
finally U is uniform on [0, 1] . 

Finally, we present the Laplace transform of . 

Proposition 6. The Laplace transform of is: 



E 



ex P (-AAy)j = (/C,(V2A)j — ^ ^ I _ 



where I u and /C„ denote the usual modified Bessel functions, with parameter v (we 
use Ylu instead of K v so that no confusion with the strike K may occur). 

3.3. Proofs of Propositions SI 13 El We follow the rationale of the proofs of 
Propositions [IJ [2j [3j after extending adequately the points (i), (ii) and (hi) in 
(|2XTjl from the BES{3) process to BES{8) process, for S > 2. 



3.3.1. Here are these extensions: 

(i),5 the Doob /i-process relationship between BES(S) and BES(4 — 5), killed 
upon hitting 0, is 



tAT ) 



1v t p(4-5). 



where 5 = 2(1 + v). 

{\i)s the particular instance of formula |T]) with M t = {X t /\T Q ) 2 ' J under P^^: 



E 



(4-5) 



(c 2 » - (x tATo f»y 



£r2i/^(4-<5) 



Ftl(r, <t) 



with 7 C = sup {t < Tq : X t — C}. 
the time reversal result: 



(X(r -t)) * — ^o) under p[ A ^ is distributed as {X tl t < gi) under Pq 



(S) 



3.3.2. The preceding results lead us to the following descriptions of r% (t) in terms 
of first and last passage times: 

rj*>(t) = P?- S \j k <t<T ) (3), 
= P^( B i>t)-Pj; Si (g 1 >t) (4), 



where k = 1/K 



1/(5-2) 



3.3.3. Proof of Proposition [5] 

It suffices to follow the arguments of the proof of Proposition replacing every- 
where dimensions 1 and 3 by dimensions (4 — S) and S. □ 
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3.3.4. Proof of Proposition [HI 
Starting from (3)^ we obtain: 



E 



(4-5) 



exp ( — \A K K 



dt e 



-A/ 



5 
PA 
5 



E (A-S) ( e -AT ) 



1 



1 



- 1 



with the help of the time reversal result (hi) s- Now, classical computations of 
Laplace transforms for first hitting times and last passage times of Bessel processes 
yield (see, e.g., Kent 0], Pitman- Yor [TO]): 



E { 5) (e-^ 1 



2A 



(V2Afc)' 



„ (V2A) 



2T(y + l)J„(fcV2A) 

where JC V and /„ denote the usual modified Bessel functions with index v. Finally, 
we have obtained the formula given in Proposition [5J 

3.3.5. Proof of Proposition [U 

For K — 0, the result follows by scaling, as in dimension 3. 
For K > 0, 

(i) We shall use formula (4)5 (to obtain the asymptotic result for r^\t), t — » 
oo) together with the formula for the distribution of a last passage time of a 
transient diffusion (see, e.g., Pitman- Yor [TU], as well as Salminen [T2l|13j ): 

s'(y)a(y) 



P x (g y e dt) 



My) 



-p t (x,y)dt. 



In our present case, we have: s(y) = — l/i/ 5 2 \ so that: 
a(y) = 1. Thus we have: 

'5-2" 



s'(y) = (S-2) , 



2// 



p). 5 \x,y)dt, 



so that, going back to the expression (4)5 for , we get: 



(*) 



5-2 



2 

Next, we shall use the following explicit formulae: (with v 

exp | | I v 



Pt(x,y) 



t o 



ds (pW(0,l) ~ P { s S) (k,l) 

ormul, 

(x 2 +y 2 ) 



(6) 



3-2 i 



2/ 



-(?) 



(see, e.g., Revuz-Yor [TT], Chapter XI). Thus: 

P W(o,i)-pW(M) 

1 / 1 \ 1 l / k 2 + i\ fk 

■ ex P ~7T - - JZ ex P o M - 



2" S "+ i r(i/ + i) ^ v 2s ; sk" * v 2 S ; " v s 
= exp (l) ( wrb) " F cxp (~fc) 7 - (l 

= exp (-i)(l)(^)(0) ^xTT^- cxp (-y^(3 

From ([6]), it now remains to study the asymptotic, as t — > oo, of: 
Making the change of variables: s = we get, from formula ([6]): 

r - (<) = Fi - exp ("2) 

Ff(iT^) _ cxp ^ (fcx) 

Letting: j/ = fcn, we get: 



k»J y'-*WfT7)-'^'^ < 7 > 



Now, we have (see, e.g., Lebedev [5]): 

1 

3=1 




"IY 

y^o \2) \T(l + v) V 2 
Finally, going back to 0, we get: 



i / k 



2 



Thus, we have obtained the asymptotic result in (i) of Proposition [4j with 
Cjp = 2 v + l (v+i)r{v) • Note that in the particular case 5 = 3, hence: v = i, 
we get: = \ —h=jfi as claimed in Proposition Q] 



(ii) Now, we study the asymptotic as t — > 0. The first result is easily obtained, 



using the convergence in L 1 of } 2 to 1, as t — ► 0. 
We now give the details in the case if = 1. As previously, we use: 



,(4-5) 



L (t<T ) 



Now, 1 — i2j?ATo> ^ ^ 0, is a martingale under p} 4 , which may be written 
as: 



(8) 



7V t M = (2^) 7/rT0ds( ^_ 1))2 
with (7 U , u > 0) a standard Brownian motion. We then write: 



ri a >(t)=25i 
We shall show 



(4-5) 



Nl 



E 



(4-5) 



A7 



(T <t) 



as well as: 

P^To < t) = P (5) ( 5 i < f) = o(Vi). 
Indeed, from ©, we deduce that the LHS of (|10p is: 



(2v)VtE 



(2i/)>/t£ [(71). 



(2i/)Vi 



2tt 



where ( 7u , u > 0) denotes a standard Brownian motion. 
Concerning (TTTj) , we use Getoor's result [3] 

(law) 

5 i = 1/27,,, 
where 7„ denotes a gamma(^) variable; thus: 



1 

i>) 



l/2t 



I Z" 00 



6 '2l +y 



(9) 



(10) 
(11) 



For v > 1, this quantity is equivalent to: r^jy 6- ^ 2 *' (^)" whereas, for 
v < 1, it is majorized by: e - ^ 1 / 24 - 1 . In any case, we easily deduce ([TT]) . 
Remark : Using (j8]) and (fTT|) . we see that: 

rW(t)~^((l-Pf) + ). 

It should be possible to use the explicit form of the semigroup density "P*(l, r), for 
5' = 4 — <5, at least when <5 < 4, to derive (fTU]) more directly. 
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3.3.6. A more direct proof of the estimate: r^p (t) = O (p^r), as t — ► oo. 
Similar to the case of dimension 3, we have, for every dimension S > 2: 

r${t)<4*((±-K) + )=4*((4--K ^ 



so that: 



Rt J J V Wi^Ri 
y/sflK 



rf{t) < c 

Jo 



< 



r \ S/tK / -, \ 



C 1 /' I 



j 2i> -r 2 /2 

ar r e ' 

2is+l 



y/t2v + l \^lK 

4. Drawing the graphs of (r$(t),t > 0) 

4.1. In order to facilitate the drawing of these graphs, we need to use the simplest 
possible integral representations of these functions. We shall rely essentially upon 
formula (JT]) which was the key of our proof of Proposition [4j 

..«>„> _ jl t" ( d JL\ .-,/» ( «" 



Using again the decomposition: 

y\*3 1 



wh„e: £«-(!)' £(§) 



V J'=1 



jir(i/ + j + i) I ' 



we obtain: 



r^(t)=rf(t)-rf(t), with 

r (<5)l^ - JL f fc /* 1 -(v/2k) ( i 

' K \ L ) — k» JO \2/ r(i/+l) V 

r^W =^/o fc/t ^*( fe ^)7, (2/ ). 



e 2 



2/ 

4.2. In the case 5 — 3, a slightly different approach leads to the following: 




and we note that: 



0{A)=N{A)-\f- (l-e^/ 2 )) 



(14) 



where: 

N(A) 

Indeed, from formula (|12[) . we deduce 
From formula (|15p . we have 



dy e 



-y 2 /2 



lM , I"* 



l-k\ l + k (l + k 



V 



2k 2 



V 



(15) 



(by (E 



(1-Jfe) (1 + fc) + 



2fc 2 



k y^tj 2fc 2 
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V 7 * 
1 - fc 



V 



i 

~2k 
N 



1 - k 



V 
V 



Vt 
1 

" 2 



9 



1 + fc 



AT 



1 - fc 



A r 



1 



(1 - fc) V 7T 

1 - k 



N 



- | | _ ( it 



2fc 2 
l + k 

~V 



1 + fc 

~V 

Vt 



V 



(16) 



V 



V 



V 



(1 - k) V 7T 



(1-fc)" 
1 — e 2t 



(1 + fc) V 7T 



v 7 / 



-|l-e-^ 



1 

2fc2 



N 
V 



1 
2 

1 - k 

~v 



N 



1 - k 

~v 



N 



(1 + fc) V 7T 
1 + fc 



1 — e 2t 



A r 



1 + fc 

"7T 



Vt 



2fc(l + fc) V 7T 



1 — e 2t 



2fc(l - fc) V 7T 
/~2 



2fc 2 (l - fc) V 7T 



1 — e 2 * 



(l-fc)^ 
1 — e 2t 



- 1 - e ~ 



2fc 2 (l + fc) V 7T 

Formula (|17| involves 5 terms with TV, and 4 terms with "exp". Thus, we write: 

V*\t)=Mt)+rt P \t), 

with 



(17) 



hit) 



N 
1 

'2k 2 



1 

TV 



1 

~ 2 
1 - fc 



iV 



l - fc 
~7T 



A r 



1 + fc 

~Vt 



-N 



l + k 

~Vt 



(18) 



;r( cx P) 



t \ 1 1 \ 1 



w - V27iU-Fj(i^) ,1 - e ~ 



1 1 \ 1 / 

1 - e at 



jfc fc 2 y (i + k) 

t ( 1 / (i-t) 2 \ 1 / (i+fcr 



We note again the further simplification of (fill)) : 

(exp),.s [% fl\ f ( l + fc 2 \\ , /'/ 



4.3. We now present the graphs of rj^-(i) for various dimensions <5 and strikes if. 
In fact, we have drawn two kinds of graphs: 

(a) Figure 1: for each dimension 6, we consider k = 1,^,^...,^, and draw all 
the graphs together. We use formula (| 16[) for 5 = 3, and formula (J7J) for 
5 = 5,7, 9, 11, 13 to draw the graphs. 

(b) Figure 2: for each k, we draw the graphs for different dimensions. Here, we 
use formula |7J) to draw graphs for all dimensions. 
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Figure 2. Varying S 
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